FGMRES for linear discrete ill-posed problems
نویسندگان
چکیده
GMRES is one of the most popular iterative methods for the solution of large linear systems of equations. However, GMRES generally does not perform well when applied to the solution of linear systems of equations that arise from the discretization of linear ill-posed problems with error-contaminated data represented by the right-hand side. Such linear systems are commonly referred to as linear discrete ill-posed problems. The FGMRES method, proposed by Saad, is a generalization of GMRES that allows larger flexibility in the choice of solution subspace than GMRES. This paper explores application of FGMRES to the solution of linear discrete ill-posed problems. Numerical examples illustrate that FGMRES with a suitably chosen solution subspace may determine approximate solutions of higher quality than commonly applied iterative methods.
منابع مشابه
Ill-Posed and Linear Inverse Problems
In this paper ill-posed linear inverse problems that arises in many applications is considered. The instability of special kind of these problems and it's relation to the kernel, is described. For finding a stable solution to these problems we need some kind of regularization that is presented. The results have been applied for a singular equation.
متن کاملProjected Tikhonov Regularization of Large-Scale Discrete Ill-Posed Problems
The solution of linear discrete ill-posed problems is very sensitive to perturbations in the data. Confidence intervals for solution coordinates provide insight into the sensitivity. This paper presents an efficient method for computing confidence intervals for large-scale linear discrete ill-posed problems. The method is based on approximating the matrix in these problems by a partial singular...
متن کاملOn the Choice of Subspace for Iterative Methods for Linear Discrete Ill-posed Problems
with a large matrix A of ill-determined rank. Thus, A has many “tiny” singular values of different orders of magnitude. In particular, A is severely ill-conditioned. Some of the singular values of A may be vanishing. We allow m ≥ n or m < n. The right-hand side vector b̃ is not required to be in the range of A. Linear systems of equations of the form (1) with a matrix of ill-determined rank are ...
متن کاملThe structure of iterative methods for symmetric linear discrete ill-posed problems
The iterative solution of large linear discrete ill-posed problems with an error contaminated data vector requires the use of specially designed methods in order to avoid severe error propagation. Range restricted minimal residual methods have been found to be well suited for the solution of many such problems. This paper discusses the structure of matrices that arise in a range restricted mini...
متن کاملFractional regularization matrices for linear discrete ill-posed problems
The numerical solution of linear discrete ill-posed problems typically requires regularization. Two of the most popular regularization methods are due to Tikhonov and Lavrentiev. These methods require the choice of a regularization matrix. Common choices include the identity matrix and finite difference approximations of a derivative operator. It is the purpose of the present paper to explore t...
متن کامل